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We analyze the statistics of daily price change of stock market in the framework of a statistical 
physics model for the collective fluctuation of stock portfolio. In this model the time series of price 
changes are coded into the sequences of up and down spins, and the Hamiltonian of the system is 
expressed by spin-spin interactions as in spin glass models of disordered magnetic systems. Through 
the analysis of Dow- Jones industrial portfolio consisting of 30 stock issues by this model, we find 
a non-equilibrium fiuctuation mode on the point slightly below the boundary between ordered and 
disordered phases. The remaining 29 modes are still in disordered phase and well described by Gibbs 
distribution. The variance of the fluctuation is outhned by the theoretical curve and peculiarly large 
in the non-equilibrium mode compared with those in the other modes remaining in ordinary phase. 



The number of days in which prices of all stock is- 
sues in Dow- Jones industrial portfolio moved in the same 
direction is 80 within 3636 trading days in the period 
from 9/Jul/86 to 22/Nov/OO, while its probability is 2^29 
when we assume Bernoulli trials. How can we explain 
the factor 10^ in the difference between these values? 
The methods and the concepts, as scaling and criticality, 
developed in material science have been applied to the 
study of financial markets |2[ || . Recently a framework 
based on spin glass model to study the collective price 
changes of stock portfolios was proposed . The ap- 
plication to 1-minute price changes in D- J portfolio made 
clear that the concept of energy works even in financial 
markets as well as the above physical concepts. Here we 
study the properties of daily price changes in D-J portfo- 
lio based on the same model and attempt an explanation 
of the factor 10^ by Gibbs factor of canonical distribu- 
tion. Through this study we clarify the applicability of 
the spin glass picture to the price fluctuations in a wide 
range of time scale and find a novel feature of stock mar- 
ket. That is a non-equilibrium fluctuation mode on the 
point close to the boundary between ordered and disor- 
dered phases. D-J portfolio has been quenched into the 
unstable region, but does not reach equilibrium. The 
variance of the fluctuation, which is physically equal to 
susceptibility and is called (the square of) risk in financial 
economy, is outlined by the theoretical curve and pecu- 
liarly large in the non-equilibrium mode compared with 
those in the other modes remaining in ordinary phase. 

In this paper, we concentrate on the statistics of the 
sign of price change |Q, 1). The time series of price 
changes are coded into the sequences of up and down 
spins Si — ±1 (i=l to portfolio size N) and the Hamilto- 
nian is expressed by long-range spin-spin interactions as 
Sherrington-Kirkpatrick model of spin glass Q , which is 
given by 



H[S, h] — — JijSiSj — hiSi, 



(1) 



<hj> 



We consider a portfolio as a subset of the whole stock 
market, and the complement of the subset works as heat 
reservoir. Various observable quantities are obtained as 



the statistical averages with Gibbs weight assigned to 
each configuration in equilibrium. The interaction coef- 
ficients Jij are constant but not fixed yet. The exter- 
nal field hi is set to be zero in the analysis of actual data 
We use Thouless- Anderson-Palmer (TAP) mean field 
theory j|, |l^ whose fundamental equation is given by 



tanh(^^ Ji^ 



J24i^-m')m). (2) 



of rrii =< Si > to determine those coefficients. The 
partial differentiation of TAP equation ^ with respect 
to hi yields the equation 



k 



A,k[m]Xkj 



(3) 



where Xij = dnii/dhj is susceptibility and 



Ay[m] = -J»j-24m,mj+%[^ Jffc(l~m|)+ _ J . 

k ^ 

(4) 

On the other hand, fluctuation-response theorem relates 
the susceptibility Xij the covariance Cij —< (Si — 
mi)(Sj — ruj) > as 



Xij — Cij. 



(5) 



Substituting the equation (|^) into (|3|), we can derive the 
relation between Jij and Cij . Interpreting Cij as the time 
average of empirical data over the observation time, Jij 
are phenomenologically determined by the equation 



Jij — 



1 + Jl-8Cr/m,r 



Amirrij 



(6) 



where C~^^ is the (i, j)-element of the inverse of the covari- 
ance matrix. In the case with w the approximation 
Jij « —C~j^ is apphcable. 

We investigate three datasets a, b and c of D-J port- 
folio consisting of N=30 stock issues in the framework 
of the model. The dataset a is the time series of stock 
prices sampled at 1-minute intervals in the period from 
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16-May-2000 to 21-Jun-2000, b is sampled at 10-minutes 
intervals in the period from l-Dec-1999 to ll-Nov-2000, 
and c is the daily price in the period from 9-Jul-86 to 22- 
Nov-2000. The time evolutions of Intel's stock price on 
each time scale are shown in Fig. |l| for an example. The 
time series of price changes are coded into the sequences 
of up and down spins. The covariance Cij for 435 pairs 
of i and j are derived from the coded data. Then the 
interaction coefficients Jij are calculated by the equation 
||l2|. The energy spectra of the system (portfolio en- 
ergy) are defined as the eigenvalues of the Hamiltonian 
H[S,Q]. The probability density of portfolio energy for 
each dataset is empirically obtained from the relative fre- 
quency during the observation times as 

AE AE 

p{E) = P{E- — <H[S,0]<E+—)/AE. (7) 

In our formalism that is given by the equation p{E) — 
n{E)e~^ / Z with the density of states n{E) and the par- 
tition function Z if the system is in equilibrium. The 
empirical probability weight p{E)n{{)) / p{Q)n{E) is plot- 
ted in Fig. H and compared with Gibbs weight as the 
theoretical prediction. In the panels a and b we see a re- 
markable fit of equilibrium line to empirical data. From 
these figures, the price changes sampled at 1 and 10- 
minutes intervals are well described by the model and are 
close to equilibrium. On the other hand, in the panel c 
some deviation of data from theoretical line is observed. 
Monte Carlo simulation of the probability distribution 
of the system magnetization ra — (l/N) X^il^i 'S'ij which 
describes the degree of the alignment of price changes in 
portfolio, make more clear this difference in the nature of 
price changes on the different time scales. The result is 
shown in Fig. |^. We see the theory explains the empirical 
data very well in the panels a and b, while the dataset c 
shows entirely different profile from the theoretical pre- 
diction. Does this difference indicate the break down of 
the theory on the time scale of day? We investigate the 
theory in more detail by the mode analysis of TAP equa- 
tion (^, and intend to explain this phenomenon in the 
framework of the theory. 

The stability of the ordinary phase is analyzed by the 
linear analysis of the equation (||). In our case, that is 

^^„[0]m, (8) 

j 

The diagonalization of symmetric matrix Aij of the equa- 
tion (|^) solve the equation (||) as m\ = h\/A\, where 
A\ are eigenvalues of Aij, m\ — ^ < X \ i > nii, and 
hx — ^ < X \ i > hi with the real orthogonal eigenvec- 
tors < A I i > . The critical temperature Tc is determined 
by the equation 

^A„„ = (9) 

A\^.^ is the minimum eigenvalue of Aij , which corre- 
spond to the maximum eigenvalue of the matrix Jy in 



the linear approximation. However f3 = 1/T is included 
into Jij in our case, so A\^^^ itself fills the role of the 
measure of the distance from the critical point. In our 
systems a, b and c, ^A„i„ = 0.353, 0.195 and -0.018 re- 
spectively. The systems c is in the position slightly below 
the critical point, while the systems a and b remain in 
ordinary phase. In order to visualize the statement, TAP 
free energy 1^ of the system c is shown in Fig. ^, in which 
the parameter space {rtii} is projected on the mode mx. 
In Fig. 1^, the upper and lower curves represent TAP free 
energy as the function of the modes m\ belonging to the 
maximum and the minimum eigenvalues of Aij respec- 
tively. The system c is in the unstable region, because 
the lower curve of the free energy is concave at rnx^.^ =- 
0.067. The behavior of the price fluctuation belonging 
to this mode is entirely different from the remaining 29 
modes. The probability of the magnetization fluctuation 
S\ ^ J2 < ^ \ > '^^ the unstable mode and that 
of a ordinary mode belonging to the maximum eigen- 
value A\^^^ are plotted in Fig. |5| by filled circles and 
filled boxes respectively. The fluctuation distribution of 
the modes in ordinary phase are narrow and shows good 
agreement with the equilibrium distribution, while that 
of unstable mode is broad and does not reach equilibrium 
distribution. 

The variance of each mode in the system a, b and c is 
plotted against eigenvalue A\ in Fig. ||. Those are com- 
pared with the theoretical curve 1/^a for infinite equi- 
librium system. The variance of the fluctuation, which is 
physically equal to susceptibility and is called (the square 
of) risk in financial economy, is outlined by the theo- 
retical curve and peculiarly large in the non-equilibrium 
mode compared with those in the other modes remaining 
in ordinary phase. Closing to the critical point A\ = 0, 
however, we find a systematic deviation of data from the 
theoretical curve, which is due to the non-equilibrium 
property and the finite size effect of the systems. 

The factor 10'' problem prompting us to introduce the 
interaction energy between stocks in the analysis of its 
collective fluctuation was partially explained by Gibbs 
weight in our model besides the single non-equilibrium 
mode. The physical picture of financial markets given 
here will be useful in the risk management and the selec- 
tion of portfolios. The applicability of spin glass model 
to economical system seems to suggest the ubiquity of 
the applications to other fields such as competitive or 
frustrated ecological systems. 
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FIG. 1; Time evolution of Intel's stock price on three time 
scales, a The data sampled at 1-minute intervals in the period 
from 16-May-2000 to 21-Jun-2000. b The data sampled at 
10-minute intervals in the period from l-Dec-1999 to 11-Nov- 
2000. c The daily data in the period from 9-Jul-1986 to 22- 
Nov-2000 (3636 days). 
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FIG. 2: Semi-log plot of the probability weight of portfolio energy E. Filled circle (•):the empirical probability weight 
p{E)n{0) /p{0)n{E) .The probability density 'p(E) is the empirical result, and the density of states n{E) is numerically ob- 



tained by 2 random sampling from 2 configurations. AE = 0.1. Solid line: Gibbs weight e 
a. b The result for the dataset b. c The result for the dataset c. 



a The result for the dataset 
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FIG. 3: Probability of the system magnetization. Filled circle (•):tlie relative frequency. Circle (o): Monte Carlo simulation of 
2^^ steps, a The result for the dataset a. b The result for the dataset b. c The result for the dataset c. 
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FIG. 4: Free energy as the function of the magnetization m\ for the system c. The upper and lower curves represent TAP free 
energy as the functions of the modes mx^ax = X) < Amax \ i > rm and m\^^^ = X) < Xmax \ i > rm respectively. 
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FIG. 5: Probability of magnetization fluctuation Sx for the dataset c. Filled boxes (■) and circles (•) represent relative frequency 
P{S - AS/2 < Sx < S + AS/2) with AS = 0.1 for the mode Sx^^^ = E < ^ma^r \ i > Si and Sx^,„ = E < Ami„ \ i > Si 
respectively. Boxes (□) and circles (o) represent Monte Carlo simulations based on the model. 
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FIG. 6: Variance of the modes Sx against the eigenvalue Ax- Circles (o), triangles (A), and boxes (□) represent the variance 
for the datasets a, b and c respectively. Solid line represents the theoretical curve 1/Ax for infinite equilibrium system. 



